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Abstract 

In supersymmetric theories the radiative corrections due to heavy states could 
leave their footprints in the flavour structure of the supersymmetry breaking masses. 
We investigate whether present and future searches for the muon and electron EDMs 
could be sensitive to the CP violation and flavour misalignment induced on slepton 
masses by the radiative corrections due to the right-handed neutrinos of the see- 
saw model and to the heavy Higgs triplets of SU{5) GUT. When this is the case, 
limits on the relevant combination of neutrino Yukawa couplings are obtained. Ex- 
plicit analytical expressions are provided which accounts for the dependencies on 
the supersymmetric mass parameters. 
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1 Introduction 



In low energy supersymmetric extensions of the Standard Model (SM), unless sparti- 
cle masses are considerably increased, present limits on flavour violating (FV) decays 
and electric dipole moments (EDMs) respectively allow for a quite small amount of 
fermion-sfermion misalignment in the flavour basis and constrain the phases in the 
diagonal elements of sfermion masses, involving the parameters fi and A, to be rather 
small. The bounds on the supersymmetric contribution to lepton (L)FV decays and 
EDMs have the advantage, as compared to the corresponding squark sector ones, of 
being not biased by the SM contribution - nor by the non-supersymmetric seesaw 
contribution j2] . Experimental limits on LEV decays and EDMs are then a direct 
probe of the flavour and CP pattern of slepton masses - see e.g. Ref. |Sj for a recent 
collection. From a theoretical perspective, understanding why CP phases and devia- 
tions from alignment are so strongly suppressed is one of the major problems of low 
energy supersymmetry, the CP and flavour problem. On the other hand, precisely 
because FV decays and EDMs provide strong constraints, they can share some light 
on the features of the (possibly) supersymmetric extension of the SM. 

Indeed, it is well known that even if these CP phases and misalignments were 
absent - or suppressed enough - from the effective broken supersymmetric theory 
defined at Mpi, they would be generated at low energy by the RGE corrections due 
to other flavour and CP violating sources already present in the theory, in particular 
from the Yukawa couplings of heavy states like the right-handed neutrinos of the 
seesaw model jl] and the Higgs triplets of SU{5) grand unified theories (GUT) 
Effects of radiative origin have the nice features of being naturally small, exactly 
calculable once specified a certain theory and - most interestingly - if experiments 
are sensitive to them, they yield limits on some combination of Yukawa couplings ^. 

Recently, it has been stressed that the present (planned) limit on /i — > 67 [H] 
(r — >■ /i7 [Zj) is sensitive to the lepton-slepton misalignments induced by the radiative 
corrections in the framework of the supersymmetric seesaw model and that the 
associated constraints on neutrino Yukawa couplings have a remarkable feedback 
on neutrino mass model building [3 . However, LEV decays cannot provide any 
information on the pattern of CP violation of slepton masses, while EDMs are 
sensitive to both LF and CP violations. Since the present sensitivities to the electron 
and muon EDMs, de < 10~^^ e cm and d^ < 10~^^ e cm |JJ^, could be lowered by 
planned experiments by up to three to five {12, 13^ and six to eight [HI El orders of 
magnitude respectively, it is natural to wonder whether lepton EDMs would explore 
the range associated to LF and CP violations of radiative origin. 

In this work we analyze the predicted range for d^ and c?^ when the radiative 
corrections due to the right-handed neutrinos of the seesaw model and, possibly, the 

^Of course, several effects of different origin could be simultaneously present, but there would 
be no reason for a destructive interference among them to occur. 
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Higgs triplets of SU (5) GUT provide the main source of CP violation and misalign- 
ment in slepton masses. This allows to extract many informations. If experimental 
searches could explore this range, one could obtain limits on the imaginary part of 
the relevant combination of neutrino Yukawa couplings. Moreover, the eventual dis- 
covery of a lepton EDM in this range might be interpreted as indirectly suggesting 
the existence of such a kind of fundamental particles which are too heavy to be more 
manifest. On the other hand, finding de and/or above this predicted range would 
prove the existence of a source of CP (and likely also LF) violation other than these 
heavy states. 

The pure seesaw case has been considered in Ref. [15] - see also the related studies 
fn\ ITH] on specific seesaw textures -, where it has been pointed out that threshold 
effects due to hierarchical right-handed neutrino masses enhance the radiatively- 
induced Im{Aii), i = e,fi,T. By extensively reappraising this framework, we find 
that for tan /5 > 10 the amplitude with a LL RR double insertion - proportional to 
tan^ j3 - dominates over the one involving Xm{Aii) - insensitive ^ to tan (3 -, the exact 
ratio depending on the particular choice of supersymmetric mass parameters. We 
provide general expressions from which it is easy to recognize the model dependencies 
and which complete previous analyses. Lepton EDMs are then strongly enhanced 
in models with large tan/3. To see how close experiments are getting to the seesaw 
induced lepton EDMs range, we compare the upper estimates for the radiatively- 
induced misalignments and CP phases with the corresponding present and planned 
experimental limits collected in 0. 

The range for a seesaw- induced d^ turns out to be quite far from present searches 
and, in particular, its eventual discovery above 10~^^ e cm would prove the existence 
of a source of CP violation other than the neutrino Yukawa couplings. On the 
contrary, the present sensitivity to de explores the seesaw-induced range in models 
with large tan (3 and small R-slepton masses, say mn around 100 — 200 GeV. The 
planned improvements for d^ would allow to test also models with niR up to the TeV 
region and moderate tan (3. When present or planned experimental limits turns out 
to overlap with these allowed ranges, bounds on the imaginary part of the relevant 
combinations of neutrino Yukawa couplings are obtained and plotted in the plane 
(Mi,m/j), respectively the bino and the average R-slepton masses at low energy. 
These plots allow to check whether any particular seesaw model is consistent with 
present data and, if so, which level of experimental sensitivity would test it. In any 
case, an experimentally interesting contribution requires hierarchical right-handed 
neutrino masses. 

This feature is no more necessary when, in addition to the seesaw, a stage of 
SU (5) grand unification is present above the gauge coupling unification scale. It is 
well known that the main drawback of minimal SU{5) is that the value of the triplet 
mass, Mt, required by gauge coupling unification is sizeably below the lower bound 
on Mt derived from proton lifetime (see for instance for recent reviews and 

^In Ref. 52] a dependence on tan /3 arises due to the particular class of seesaw textures studied. 
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references). In our analysis we therefore keep the triplet mass as a free parameter. 
We nevertheless exploit the minimal SU (5) relations between the doublet and triplet 
Yukawa couplings since in general they are mildly broken in non-minimal versions 
of SU{5). The simultaneous presence of right-handed neutrinos and heavy triplets 
turns out to further enhance the amplitude with the LL RR double insertion over 
the one with Xm{Aii), essentially unaffected by triplets. This cannot be derived by 
the - otherwise elegant - technique based on the allowed invariants ['2U\ . 

When right-handed neutrinos and Higgs triplets are simultaneously present, the 
predicted range for the radiatively-induced de is already sizeably excluded by the 
present experimental limit. This in turn is translated into a strong constraint on 
the imaginary part of the combination of Yukawa couplings which is relevant for the 
LL RR amplitude. Such a constraint could be hardly evaded even for large values of 
the triplet mass and unfavorable supersymmetric mass parameters. The radiatively- 
induced in the presence of triplets should not exceed 10~^^ e cm, as was the case 
for the pure seesaw. However, at difference of the latter case, in the former one o?^ 
only mildly depends on the spectrum of right-handed neutrinos. Notice that planned 
searches for d^ could constrain (depending on the triplet mass mass, of course) the 
imaginary part of the combination of Yukawa couplings whose absolute value could 
be independently constrained by the LFV decay r — /X7. 

The paper is organized as follows. In Section 2 we introduce our notations, dis- 
cuss the framework and make some preliminary considerations. Section 3 considers 
the pure seesaw case by separately analyzing the flavour conserving (FC) and flavour 
violating (FV) amplitudes contributing to lepton EDMs. In Section 4 the frame- 
work of seesaw and SU (5) is discussed along the same lines. Concluding remarks 
are drawn in Section 5. Finally, in Appendix A and B we collect the RGE in the 
case of the seesaw, without and with a minimal SU{5) unification respectively. 

2 Framework and Method 

In this section we recall the expression for the lepton EDMs in the mass insertion 
approximation 1221 12ni 12^1 0] to display the supersymmetric mass parameters 
that are constrained by the present and projected searches for de and d^. We then 
draw some preliminary considerations to introduce our procedure to calculate the 
radiative 1-loop contribution to these mass parameters. The relevant RGE can be 
found in the Appendix. 

We adopt here the following conventions for the 6x6 slepton mass matrix in the 
lepton flavour (LF) basis where the charged lepton mass matrix, m^, is diagonal: 



where is the 3x3 matrix of the trilinear coupling, the A— term. All deviations 




(1) 
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from alignment in this mass matrix are gathered in the 6 matrices, which contain 
30 real parameters (including 12 phases) and are defined as: 

A\vd — /i tan (3mi = {Al^Vd — /i tan (3mi) + 1711171^6^^ (2) 

where mL^m^ are average masses for L and R sleptons respectively and ~ 
Oinisusy'fnfjvd) are the diagonal elements of Ag, so that 5^^ has only non-diagonal, 
fiavour violating, elements. 

The supersymmetric contributions to dj, (z = e,yU, r), can be splitted in two 
parts, involving respectively only fiavour conserving (FC) or fiavour violating (FV) 
elements of the slepton mass matrix (|T|): 

= df^ + df (3) 



= L I ?2 \n [mam{^,)t^n(3{lB + -lL-lR + l2)-vdIm{Al)lBm 
2 4-7r|/ip cos^ t^VK 2 

dr = L^l^p^s^g^ K^.(Jm(^^^^^^),4,, +Tm(^^^^^^),/l,^^) 

+ tan/?(Xm(5^^/ir/,m,5^^),,/^ + Jm(5^^/i*r7>,5^^),,J^^)] (5) 

where rji = i3 — A^Vd/ {mifita.n j3) , the functions / are defined as in pj and terms 
that are less relevant ^ or higher order in the 5's matrix elements are omitted. 
Notice that rje ^ i3 for relatively large values of tan/3, favored in mSUGRA and 
for which the LF and CP violations are most likely to be detected. The FC and 
FV contributions could result from different seeds of CP violation but could also 
be correlated in many different ways in models. Anyway, the experimental limit 
can be put on both because, due to the different nature of the many parameters 
involved, an eventual cancellation between these contributions appears unnatural. 
The distinction between the FC and FV contributions is also phenomenologically 
relevant: some of the |(5|'s in the FV terms are already constrained to be smaller 
than 0(1) by LFV decays, while Xm(/i) and Xm{Aii) are directly constrained by the 
EDMs. 

Our aim here is to estimate the radiative contribution to the lepton EDMs in- 
duced by the seesaw interactions, first alone and subsequently accompanied by a 
stage of SU{5) GUT. Since the present experimental bounds on LFV decays and 
EDMs already point towards family blind soft terms (i.e. sparticles with the same 
quantum numbers must have the same soft terms) with very small CP phases, at 
the cut-off A = Mpi corresponding to the decoupling of gravitational interactions 
we assume real and flavour blind soft terms, namely in eq. (P), 

^II = ^RR = "^oi3 , Ae = y^ao , (6) 

•^E.g. a contribution to df^ of the form of with 5^^ — > {A'^^Vd — ^ita,n f3mi)/ {mLmn), is in 
principle present, but it is neghgible with respect to Q). 
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with real Oq, mo and /x term. In the next sections we separately study the FC 
and FV contributions to di and obtain explicit approximate expressions for Xm{Aii) 
and the products of (5's in (0). The effects of a more general family independence 
assumption are important but not crucial and can be easily included in our analysis. 
By means of these general approximated expressions we will: 

i) derive the upper prediction for the radiatively induced leptonic EDMs, stressing 
the model dependences; 

ii) compare it with the experimental limits; 

iii) when allowed by the experiment, obtain an upper bound on the imaginary part 
of the relevant combination of Yukawa couplings. 

A couple of preliminary considerations are in order before presenting the results 
of the next sections. 

2.1 On the naive scaling relation 

In the limit that all slepton masses are family independent, the FV contribution 
vanishes and the FC one is proportional to the mass of the z— th lepton {Xm{Aii)vd ~ 
Xm(ao)mj) leading, except an accidental cancellation ^ with the yU-term amplitude, 
to the "naive" scaling relation 

di/dj = mi/rrij . (7) 

Then, due to the present experimental limit on dg, d^ could not exceed dem^/rrie ~ 
2 ■ 10~^^ e cm, which roughly corresponds to the planned sensitivity and, if the limit 
on de were still to be lowered, next generation experiments would have no chance of 
measuring rf^. 

Such considerations could provide interesting informations because (jT)) strictly 
apply only to the FC yU-contribution to the EDM, while in general both Xm{Aii) and 
the FV terms may strongly violate it. This is the case for the radiatively induced 
1m{Aii). Some of the FV contributions are instead naturally proportional to a 
different lepton mass, m^, possibly heavier than rrij, as discussed in - and, 

before, for the quark sector, in j2ni- In particular, it will turn out in the next sections 
that the FV contribution can even take over the FC one. Hence, a value of d^ above 
~ 2 ■ 10~^^ e cm is a possibility that deserves experimental tests and, interestingly 
enough, it would imply the source of lepton EDM being either the FV contribution 
or a non-universal Xm(y4jj), so providing a remarkable hint for our understanding 
of CP violation. 

''For dedicated studies on cancellations between amplitudes, also in more general frameworks, 
see e.g. PSIEI) . 
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2.2 On the relevant combination of Yukawa couplings 

Once a theory is specified, it is relatively easy to list the combinations of Yukawa 
couplings appearing in the slepton mass radiative corrections that may contribute 
to lepton EDMs. However, the actual calculation is more involved as we now turn 
to discuss. 

As an example, let us look for Im{Aii) by studying the evolution of Ae in the case 
of degenerate right-handed neutrino masses, M, for simplicity. Let us first consider 
the case of the pure seesaw and define ^ the hermitian matrices E = ylue, N = yluu- 
Notice that E is real and diagonal in the LF defining basis and N is diagonalized 
by a unitary matrix similar to the CKM one with only one phase (even for non- 
degenerate right-handed neutrinos). By solving the RGE for A^., eq. (j37|) . linearly 
in ^3 = l/(47r)^ ln(Mp;/M), Af, can only be proportional to y^E and yeN, whose 
diagonal elements are real. At 0(^3) only yeEE,yeNN,yeNE and yeEN appear, 
whose diagonal elements are again real. A potential Im^Au) shows up at 0(^3), 
through Im{yeN[N, E]N)ii. On the other hand, when also Yukawa interaction of 
the SU{5) triplets are present, eq. (jlH|) shows that, at first order in It = l/(47r)^ 
ln(Mp;/Mr), Ae can be proportional to yeE, yeN but also to U*ye, where U = ylyu- 
The latter have real diagonal elements but allow at 0(t|.) for a combination with 
diagonal imaginary part, namely {U*yeN)ii. 

Accordingly, in order to evaluate the actual coefficients in front of the prod- 
ucts of Yukawa coupling matrices which are likely to have phases, we solve the 
RGE for the soft parameters by a Taylor expansion in the small parameters tj/ = 
l/(47r)^ ln((5i/Q/) associated to the intervals between the successive decoupling 
thresholds of the various heavy states. For instance, by integrating the RGE for 
Ae in the case of SU (5) plus seesaw, a non vanishing coefficient is obtained for the 
combination Xm{U*yeN)ii at 0(t|.). Now, lepton FV transitions and CP phases 
are naturally defined in the LF basis where ye and the Majorana masses are 
diagonal and real. This basis is not invariant under the RGE evolution and one 
should diagonalize ye and Mr again at the lower scale. Therefore, one has to find 
out the effect of these final rotations. 

We adopt the rotating basis method introduced in Ref. |2S1, where the RGE are 
modified to incorporate the fact that the matrices are defined in the LF basis at 
each scale, so that ye and Mp are always diagonal. It is worth to stress that, even 
if the rotations needed to diagonalize ye are small, their effect could be crucial for 
Xm{Aii) and the four products of 5's in eq. (0). For instance, in the case of SU{b) 
plus seesaw it turns out that this correction exactly cancels the term proportional to 
t^Im{U*yeN)ii in Ae- On the contrary, these rotations can be safely neglected when 
deriving approximate expressions for the radiatively-induced LEV decays because 
the latter are only sensitive to absolute values of (5's. 

^Here and in the following Dirac mass terms are always written as fumffL- 
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3 Lepton EDMs and Seesaw 



In this section we consider the predictions for lepton EDMs in the context of the 
super symmetric extension of the seesaw model, namely we assume the MSSM sup- 
plemented with the seesaw Yukawa interactions and Majorana masses for the right- 
handed neutrinos as the effective theory valid up to the cut-off A = Mpi where 
gravitational interactions decouple. Starting with real and universal boundary con- 
ditions at Mp/, we solve the RGE displayed in Appendix A by expanding in the 
small parameters defined by the right-handed neutrino thresholds 

1 A 1 , M3 I M2 

^^-(i^O^^^Ms ^^-(i^o^^^M^ ^'-Jm^^^'W, ' 

with the ordering M3 > M2 > Mi. We thus obtain approximate analytic expres- 
sions for the radiatively induced (5's and Xm^Au) which depend on mo, ao and the 
Yukawa couplings defined in the LF basis at the scale A. Of course the latter can 
be immediately translated into the corresponding ones defined at any scale, e.g. Mi 

or fTlgy^gy. 

Precisely because of the lepton flavour and CP violating Yukawa couplings, the 
LF basis is continuously rotated and rephased with the RGE evolution and, as 
already discussed, we handle this by working in a rotating basis. It turns out that 
the basis transformation introduces negligible corrections in the FC terms (@]) but 
important ones in the FV terms Notice also that the seesaw effects stop at 
the decoupling threshold of the lightest right-handed neutrino. Mi, and that the 
RGE evolution of these effects down to the supersymmetric scales where CP and 
lepton FV transitions are estimated is generically small and can be neglect for our 
estimates. We now study separately the FC and the FV contributions. 



3.1 Flavour conserving contribution 

Starting in the LF basis at Mpi from A^, = yeCio, the seesaw interactions generate a 
Xm{Aii) in the LF basis at Mi. At leading order in the fs and defining ylPaHu = Na 
(a = 1,2,3) and the right-handed neutrinos projectors Pi = diag(l,0,0), P2 = 
diag(l, 1,0), P3 = i3, the latter reads ^: 

Im{Au) = Saoz/e. {t2t3lm{N2Ns)u + titsIm{NiNs)ii + tit2lm{NiN2)ii) , (9) 

where the Yukawas in the r.h.s. are evaluated at A. A similar ^ formula were previ- 
ously presented in Refs. (THl [13 HH] , but with the various Yukawas involved defined 

^Actually, Q is an approximation by excess and improves as the involved Yukawa couplings 
become small. However, for our estimates it is reliable up to 0(1) Yukawa couplings. 

^As far as the comparison with the expressions in Refs. ^lEIEl is possible, we find agreement 
with the last papers up to the coefficient of Xm{NiN2)ii- 
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at different scales. Then, one must be more cautious in drawing general conclusions 
and the authors validate theirs with specific numerical examples. Our approach 
offers the advantage that the corresponding results become more transparent and 
allow for an easier estimate of the effects once a pattern of seesaw parameters is 
assigned. Anyway, the crucial point fH! is that the more right-handed neutrinos are 
hierarchical, the more the FC contribution increases. Indeed, it vanishes in the limit 
that right-handed neutrinos are degenerate, in which case a contribution to Xm{Aii) 
only appears at fourth order, proportional to 2m{yeN[N, E]N)ii. Notice also that 
the naive scaling relation is generally violated jTHj by (jHl). 

Eq. Q displays a linear dependence on the unknown parameter oq. By defining 
An = |A,,|e*«<^-^s then 0a, ~ t2t3lm{N2N3)ii+ tihXm{NiN3)u+ tit2lm{NiN2)u is 
completely specified by the seesaw parameters. However, as appears from eq. 
the experimental limit on di doesn't probe directly 0^^, rather it gives a bound on 
1m{Aii)vd = miXm{ai) once supersymmetric masses are fixed (and up to unnatural 
conspiracies between the various amplitudes). Figs. [TJ taken from Ref. [S], show 
the present upper limits on \Im{ae)\/mii and the planned ones for \Im{a^)\/mR in 
the plane (Mi,m/j), respectively the bino and R slepton mass at nisusy Notice that 
these limits are quite model independent because, apart from ttir and Mi, there is 
only a mild dependence on m^, which we have fixed as in mSUGRA for definiteness. 
Indeed, since the A-term amplitude in df" arises from pure bino exchange, it does 
not involve * /i nor tan/3. In mSUGRA there is an unphysical region in the plane 
(Ml, m/j) corresponding to rn^ < and which has been indicated in light grey in the 
plots. Anyway, the dark grey region and below is also excluded because < Mi, 
in contrast with the requirement of neutrality for the LSP. 

To find an upper estimate for the seesaw induced df" , let us evaluate \Tm{ai)\/mR 
from eq. by considering only its first term. This situation is representative be- 
cause the terms proportional to ti are negligible when the lightest right-handed 
neutrino has smaller Yukawa couplings, as happens in many models and as one 
would guess from similarity with the charged fermion sectors. Anyway, the follow- 
ing discussion is trivially adapted to other cases. For definiteness, we adopt this set 
of reference threshold values: Ms = 10^^ Q^y, M3 = 10^^ GeV, A = 2 lO^^ GeV. 
Since, as demanded by perturbativity, {N2N3)ii < 0(1), 

0(0.02^)^ . (10) 
rriR V 2 10 V 

This upper estimate is easily adapted to any given model once the relation between 
Oq and rriji is made explicit. To make the dependence more manifest, let us introduce 
the following two mSUGRA situations: a) Oq = 2mQ; b) = My^., where M1/2 is 
the universal gaugino mass at the gauge coupling unification scale. In the first case 



^The function Ib in eq. contains a factor that cancels the one in the overaU coefficient. 
The exact expressions for Is can be found in |3j, as well as approximate ones suitable for various 
pattern of supersymmetric masses. 
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Experiment 
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Figure 1: Experimental upper bound on \Tm{ae)\lmR and \Tm{a^)\/mR for de < 
10"^'' e cm and < 10~^^ e cm respectively jB]. corresponds to the slope: 
r = Mi/rriR. For the present sensitivity d^ < 10~^® e cm, the numbers have to be 
multiplied by 10^. mj;, is fixed as in mSUGRA. 



{ao/mji)^ ^ 2(1 — 0.9r^), with r = Mi/rriR. The ratio ao/niR thus displays a mild 
excursion as it decreases from \/2 down to 0.45 when moving from the vertical left 
axes, where r -C 1, to the joining line of the dark grey region, where r = 1. For case 
b), the situation is opposite and not so mild: {ao/rriji) ^ 2.5 r. In the more natural 
cases in between, the ratio ao/m^ should thus be rather stable. 

Let us now compare the upper estimate ()10|) with the experimental bound. As 
appears from eq. the latter improves linearly with the experimental sensitivity 
to di. Taking ao ~ ruR, figs. ((H) show that values around 0.02 for \Tm{ai)\/mR 
would require an experimental sensitivity to de and d^ respectively at the level of 
X0~28 _ xo~^^ e cm and 2 10~^^ — 2 10~^^ e cm, the exact value depending on 
the particular point of the (Mi,m/{) plane. However, to avoid charge and color 
breaking, the constraint ao/mR < 3 in general applies. Thus, focusing for instance 
around niR ^ 500 GeV, the FC contribution cannot exceed ~ 10~^^ e cm for de and 
~ 2 10~^^ e cm for c?^, even with highly hierarchical right-handed neutrinos. Notice 
that the former value could be at hand of future experimental searches for de, while 
the latter is at the very limit of the planned experimental sensitivity to c?^. Allowing 
for smaller ttir values, tjir ^ 200 GeV, the FC seesaw upper estimate comes close 
to the present limit for de while it cannot be more than ~ 2 10~^^ e cm for d^. 
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The result for (i^ is thus a kind of "negative" one, since its eventual future dis- 
covery above ~ 2 10~^^ e cm couldn't be attributed to the radiative FC contribution 
of the seesaw. Even if the situation for df" appears more optimistic, it is worth 
to underline that the previous upper estimate actually applies to models with at 
least four 0(1) neutrino Yukawa couplings and large CP phases, like those in Ref. 
[IH]. However, as we now turn to discuss, the FV contribution, underestimated by 
previous analyses, could drastically enhance the di upper estimate. 



3.2 Flavour violating contributions 

In the following we isolate the potentially most important products of 5's and give 
an estimate of their relative magnitude with respect to the FC amplitude. The 
relevant approximations for the flavour violating elements, i j, of the 6's are: 



rriLmRSf^^ = aom,[-2t3Afij + ^^FA{a,a)ij + ^tatbFA{a,b)ij] (11) 

a a>b 
a a>b 

mlS^"- = -i6ml + 2al)hMj + J2^FLia,a)^j + J2iatbFLia,b)ij (13) 



a>b 



where a, 6 = 1, 2, 3 , the matrix Af is defined as tsAf = t^N^ + ^2-^2 + ^1-^1 and 

FA{a,b) = 15{E,Na}-5{E,N,} + 12{Na,N,} + A[N,,u'^^^] 

+ 2((iV,(,) + D,)Na + (iV,(,) + D,)m) + 4(D(")/ao + Di-^)N, 

+ [Na,E]+A[N,,Na] + 7[E,N,] , (14) 

Fnia,b) = 8[iQml + Aal)y,Nayl-iQml + 2al)y,N,yl] (15) 

FL{a,b) = 2{6ml + 2al){{3Nt + E,Na} + 2Di''^N,+ [Na,u'^E^]) + 2m^'j,^N, 

+ 4al{{3N, + E,Na} + {E,N,}) + 2{GL + AaoDi'^^)N, (16) 

with the Yukawas in the r.h.s. evaluated at the scale A. The subscript (d) indicates 
to take only the diagonal elements of the matrix, is defined through [u^^\ E] = 
{E, 3E + Na + Df.} and the definition of all the other quantities can be found in 
Appendix A. 

By means of the above expressions, the predictions for the imaginary part of the 
various product of 5's can be studied. An imaginary part in the products {S^^S^^)ii, 
{6^^6^^)ii, only arises at third order in at least two different t's. 
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3.2.1 The LL RR contribution 



Let us firstly discuss the LL RR double insertion, Tm{6^^ ^riem£6^^)ii, which turns 
out to be the quantitatively most interesting one. Since in general rje ^ i3 and the 
phase of /i - if any - is extremely small so that it can be safely neglected in the 
discussion, Im{5^^ iir]im(,5^^)ii ~ film{6^^m£6^^)ii. From eqs. (fT^ IT^ it appears 
that 6jij^, Sjlj^ are respectively of second and first order in the t's. Thus, the lowest 
order is the cubic and the LL RR contribution reads: 

Jm(5^^m,5««)., = g^ /^^o + 2ag)(6mg + 3ag) ^^^^^ ^ U)Xm{N^EN,)u . 



a>b 



(17) 

Notice that 2m{NaENb)ii = tan^ f3X'm{NamjNb)ii/rn'f, where rrit is the top mass. 
Due to the hierarchy in mg, a potentially important effect could only come from 
Im{NaisNf,yi) . Notice also that, on the contrary of the FC contribution, (jl7j) doesn't 
strongly depend on Oq. 

The relative importance between this amplitude and the FC one can be easily ap- 
preciated by neglecting the terms proportional to ti, which, as already mentioned, is 
justified when the lightest right-handed neutrino has the smallest Yukawa couplings: 

^FV.LRu ^ ^tan^/? (6mg + 2ag)(6mg + 3ag) ImiNsmjN^)^ I'^ 

For realistic values, the ratio of the two loop functions is slightly smaller than one. To 
obtain a rule of thumb, let us take mo ~ Oq ~ ""^l,/? and the reference ratio A/M2 = 
2 10^. Then, unless ad hoc fine-tunings in the structure of A^, 1m{N^m\N2)ii ~ 
m'^Xm{N^N2)ii-, and one finds 

- 0.5 -j^ — -^ (rule of thumb), (19) 

where |/ie«,P ~ 0.5m|. + 20M^ is the value of /i accounting for radiative electroweak 
breaking in mSUGRA. Despite being of third order in the t's, the FV amplitude can 
take over the FC one thanks to its tan^ (3 dependence. The precise value of this ratio 
is displayed in figs. |21in the plane (Mi, m/j). Cases a) and b) are separately displayed 
so that the behavior for any situation in between can be easily extrapolated. The 
relevant generalizations are also reminded. The plots show that the rule of thumb 
is quite reliable and allow to extract, for each point of the plane, the value of tan/3 
for which the FV amplitude takes over the FC one. In both cases a) and b), this 
happens for tan/? ^ 10 - the only exception being the region with r 1 for case a), 
where tan /3 > 20 is required. 
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Theory: -^^-tf^ — foi" Seesaw 




Figure 2: Ratio df^^^^^/df^ when Im{N3mjN2)ii ~ ni^Tni{N^N2)ii- We have 
assumed as reference: A/M2 = 2 10^, rriL as in mSUGRA, = /ig^ and tan/5 = 10. 



3.2.2 The other contributions 

The amphtudes with the other products of 5's are less important than the LL RR 
one. Consider for instance Xm{6^^6^^)ii. The flavour violating elements in 6^^ and 
are respectively of second and first order, so that the product appears at third 
order: 

Tm(5^^5^^),, = 8 ^-^°(l^;"o + 6ag) y. ^^^^^^^ ^ t,)lm{N^EN,)u . (20) 

a>b 

It vanishes in the limit ao — > 0, as the FC contribution. Neglecting the terms 
proportional to ti, the ratio of the LR RR amplitude and the FC one reads 

It is easy to check that this ratio is smaller than one (unless ad hoc fine-tunings 
in the structure of N): the ratio of the two loop functions is slightly smaller than 
one for realistic values of the supersymmetric parameters and, taking for instance 
mo = ao and A/M2 = 2 10^ (jHI) should not exceed ~ lO^^tan^/?. 
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For Xm{6^^6^^)ii, no imaginary part can arise at second order in t's because both 
6^^ and 6^^ are proportional to t^jV. At third order there are many contributions 
and it is lengthy but straightforward to check that they are proportional to at 
least two different fs. This contribution is also proportional to oq and could be 
comparable to the FC one but is in any case smaller than the LL RR one. The 
expression for the double LR insertion is also quite involved. It is proportional to 
Oq and, being also suppressed by a factor m^/m| with respect to the LL RR one, it 
can be safely neglected. 

3.3 Predicted range and constrains on Yukawas 

For values of tan/3 > 10, for which the EDMs are enhanced and thus most likely 
to be observed, the FV amplitude with the LL RR double insertion is generically 
dominant with respect to all other amplitudes. Then, when tan/3 > 10, all the 
considerations made previously for df-^ actually apply to di when strengthened by 
the rule of thumb factor (fT^ . To give an example, if ^2 + ^3 ~ 0.1 and fj, ~ fj,ew, 
4 cannot exceed ~ 0.5 10-28(-27) tan^/^/lO^ e cm and d^ ~ io-26{-25) tan^/^/lO^ 
e cm when rriR ~ 500(200) GeV. Planned experimental sensitivities to c?^ could 
then test the seesaw radiative contribution for models with small m^j and/or large 
tan P, in which case the factor (fT9|l could be up to ~ 50 — 60 so that dfj, should not 
exceed O(10~^^) e cm. On the contrary, the range of the seesaw induced de already 
overlaps with the present experimental limit for values of m^j up to 1/2 TeV when 
tan P ~ 30. 

Barring unnatural cancellations, planned searches for de could thus test each 
term of the sum in (|17|) , namely each ^ 

2 

tatb{ta + tb)Im{Na^Nb)u (a > 6) . (22) 

The effect of an eventual two orders of magnitude improvement for dg on the upper 
limit on Xm{Nam'jNh)ii/ml is displayed in fig. El by taking as reference values 
tan/3 = 30 and tatbita + tb) = 2 10""^. It turns out that planned limits could be 
severe enough to test models with hierarchical neutrino Yukawa couplings. This 
cannot be done by present limits. For any given seesaw model, it is straightforward 
to extrapolate from the plot the level of experimental sensitivity required to test it. 
Notice also that for d^, limits on Im^Namj N^) 22 strong as the present ones 
on Tm{NamjNh)u/ml would require a sensitivity to d^ at the level of 2 10~^^ e cm. 

^Needless to say, this remains true even if the FC were the dominant amphtude, in which case 
the bound would be stronger. 
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Upper bound on Xm{Na^Nb)ii 




100 200 300 400 500 
Ml (GeV) 



Figure 3: Upper bound on Tm{Nam\Nh)ii/m'^. Solid (dashed) lines refer to case a) 
(b)). 

4 Lepton EDMs, Seesaw and SU{b) Triplets 

We now add to the supersymmetric seesaw model a stage of a minimal SU (5) GUT 
above the gauge couplings unification scale, Mqut ~ 2 10^^ GeV. Namely, we include 
the contribution of the higgs triplets Yukawa interactions to the RGE evolution of 
slepton masses from Mpi down to their threshold decoupling scale Mr, which is very 
likely to be bigger than M3. Notations are defined in Appendix B. 

It is not restrictive to work (at any scale) in the basis where (j/J =)ye and Mr 
are real and diagonal and Uu = V^du4>uV, where du are the (real and positive) 
eigenvalues of V is the CKM matrix in the standard parameterization (more on 
this later) and (pu is a diagonal SU{3) matrix. The RGE for the radiatively induced 
misalignments are written in eqs. ()47|) to ()50p. At first order in ^1,2,3 defined in eq. 
dHl) and 

their solutions at the scale Mi reads: 

ml6^^ = -{6ml + 2al){2Uyl + 3tTU*) 

mi5^^ = -{6ml + 2al){{3tT + h)yl + tsAf) (24) 
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where the matrix Af is defined as t^N = t^N^ + ^2-^2 + ^1^1 • It is understood 
that all the quantities in the r.h.s. of are evaluated at A. The small effects of 
the subsequent evolution from Mi down to mgusy can be neglected in the following 
discussion. We explicitly write only the first order terms in the t's because, 
contrarily to the seesaw case, they already produce a potential imaginary part for 
the FV contribution, eq. Instead, for the diagonal part of Ae this is not the 

case, as we now turn to discuss. 

4.1 Flavour conserving contribution 

As anticipated in the simplified discussion of Section |21 a potential candidate for 
Xm(y4ji), proportional to Xm{U*yeN)ii, could show up at Oit"^). It is lengthy but 
straightforward to see that such term is exactly canceled by the effect of rotating 
the basis. For the same reason, in the general case with different thresholds Mi^2,3, 
Mt, the overall coefficient of Xm{U*yeAf)ii is zero. Therefore, the second order 
contribution to Jm^Aa) is just 

ImiAu) = 8ao Ve, [ + ^3)^2 ImiN^Ns)^ + tr + ^3)^1 Im{N^N3),i 

+t2tiIm{NiN2h ] , (25) 

namely the pure seesaw one discussed in the previous section, eq. (fTUI) . with the 
substitution (3/2 tx + t^), due to the fact that above also the triplets 

circulate in the loop renormalizing the wave functions, is naturally expected to 
be small (triplets will not decouple much below 10^^ GeV) so that eventual higher 
order contributions involving It are expected to be negligible with respect to (j^ . 
As a result, a stage of S'f/(5)-like grand unification, cannot enhance by much the 
FC contribution with respect to the pure seesaw case. 

4.2 Flavour violating contribution 

On the contrary, products of two 5's have an imaginary part proportional to t^tT 
Tm{U*yeN')ii and the FV contribution to di is potentially bigger than the FC one. 
Most interestingly, the predicted range for the radiatively induced turns out to 
have been already sizeably excluded by the present experimental bounds. Planned 
searches for d^ would also get close to test the range corresponding to radiatively 
induced misalignments. 

Let us consider in turn the predictions for the imaginary part of the products of 
(5's, eq. (0), from their expressions given in eq. (j24j) . As before, the most important 
contribution comes out from the LL RR double insertion. Since rji ^ i3 and the 

^°0f course eqs. H24II overestimate the misalignment. However, for our estimates here they are 
reUable up to yt ~ y,^^ ^ 1- 
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phase of /X - if any - is experimentally small enough to be safely neglected in the 
present discussion, one has at the lowest relevant order in the fs 

Xm{5'^''m,5''%i = 3 trh ^^""o + ^^"^^ Xm{MmtU*)ii . (26) 

mym'i 

This FV contribution is potentially much bigger than the FC one because, apart from 
Yukawas and numerical coefficients, it is enhanced by a factor (m^/i tan /5)/(mjao)- 

The other FV combinations in (0) are: 

= 12 t,t, 4 ^'"(^f ^•)- (27) 



2 

m. 



TTl 

Im{6^^6^^)ii = -^Im{6^^6^^)ii . (29) 



The only contribution which does not vanish in the limit ^ is the LL RR 
one. To compare the amplitudes, (j^ . (jTTj) and (j^ . have to be multiplied 
respectively by /itan/5 and rriLmji and also by the appropriate loop functions, which 
have the same sign and in general are of the same order of magnitude (for more 
details see |3J). Then, if > the four FV amplitudes have the same sign. However, 
the double LR insertion is always negligible with respect to the LL RR one because 
of the suppression factor mj/mj^. For the other amplitudes (j^ . lj^ the suppression 
factor with respect to is mij^i/(/itan/?) (actually smaller due to the numerical 
coefficients). Then, even in the case of < 0, a reduction of the LL RR amplitude 
due to accidental cancellations seems unrealistic. Of course, also the contributions 
due to different thresholds in the right-handed neutrino spectrum are present, in 
exact analogy to what has been discussed in the previous section. 

It is instructive to focus on the magnitude and dependencies of the combination 
Jm{AfmiU*)ii. For and d^, neglecting subleading terms proportional to y^, y"^ 
and defining Vtd = |Vtd|e"^*'^: 



(30) 



m 



Im{ArmeU*)n ^ yt\Vtd\ Jm(e-^<^"M3) + — Vt, Im{e-'^^^ (31) 

\ J 

where we exploited the fact that Vts is real in the standard parameterization. The 
latter is convenient to stress that the CP phases involved in the above combinations 
could be naturally large - as is indeed the case for (pt^ - but, of course, any other 
choice must give equivalent results The contribution proportional to Xm(A/i2) 

^^Had we exploited the freedom of parameterizing V in such a way that Vu are real numbers, 
then the CKM phase of V would have been hidden in the redefinition of N . 
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has important suppression factors. Moreover {N'ul is independently constrained to 
be quite small from the present limits on — > ej. A plot of the present upper limit 
on IC12I, with C = {inYt^Af, in the plane (Mi, rriR) can be found in Ref. p. The 
limit were derived for the seesaw but also applies without significant modifications 
to the case of SU{5) plus seesaw. As a result, once fixed My, experimental searches 
for de and c?^ represent a test for Xm(e~*'^" A/13) andXm(A/23) respectively. Although 
present searches for t ^ n'j {t ^ 67) are not able by now to interestingly constrain 
1^*231 (IC'isI), eventual experimental improvements would have an impact on (de) 

m 

Notice also that the naive scaling relation is violated according to 

4 ^ \Vtd\Tm{e-'^^'^Ari3) 

d, Vts Im{Af2s) ^ ^ 

and that the combinations of Yukawas relevant for di are independent on the phases 
of the diagonal SU{3) matrix, On the contrary, the latter affects (see for instance 
Ref. j22j) the proton decay lifetime due to d=5 operators, whose most important 
decay mode in the case of minimal SU{5) is p ^ K^v. 



4.3 Predicted range and constraints on Yukawas 

To understand how close to the experimental sensitivity is the radiatively induced 
EDM range, we plot in figs. El and |H1 the upper estimate for the most important 
products of (5's. We consider a degenerate spectrum of right-handed neutrinos to 
pick out just the effect of the triplets, substitute the Mp;- values for yt (~ .7) and the 
relevant CKM elements and choose as reference values K/Mt = 10^ and A/M3 = 
2 10^. Then, the upper estimate follows by requiring perturbativity, Xm{N2z) < 1, 
Xm(e~*''^"A'i3) < 1. Solid and dashed lines refers to cases a) and b) respectively. We 
don't show case b) for Im{S^^miS^^)ii/mr, because the predicted value is essentially 
fiat. The upper estimate for Xm{6^^6^^)ii is not shown, being closely related to 
that on Im{6^^S^^)ii (see eq. (f^^ ). The bounds on Im{5^^mi5^^)ii/ mr are not 
displayed because they are too small to be of potential interest. 

For an easy comparison, we have taken from the sleptonarium [3^ the exper- 
imental limits on the same quantities, figs. |3] and [7| . The experimental limits 
on Xm{6^^6^^)ii are close to those on Im{6^^6^^)ii and mildly depend on m^, 
which has been fixed as in mSUGRA in the plots. On the contrary, the limits on 
Xm{5^^mi5^^)ii/mr are proportional to l/(/itan/?). For definiteness, tan/5 = 10 
and /i = yUeto have been assumed. The experimental bounds shown in figs. 0] and 
|71 correspond to the present bound de < 10"^'' e cm and to the planned sensitivity 
d^ < 10"^'* e cm. Since the experimental bounds are proportional to the bound on 
di, it is straightforward to extrapolate the sensitivity to and de required to test 
the radiatively induced de and (i^. 
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Experiment 




100 200 300 400 500 100 200 300 400 500 



Ml (GeV) Ml {GeV) 

Figure 4: Experimental upper bounds j3] on various products of 5's corresponding 
to the present sensitivity d^. < 10~^^ e cm. 



Theory: SU (5) with Seesaw 
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Figure 5: Upper estimate for various products of 5's in SU{b) with seesaw. The 
reference values K/Mt = 10^ and A/M3 = 2 ■ 10'^ have been taken. Solid (dashed) 
lines correspond to case a) (b)). 
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Upper limit on 

Xm(e-*'^*'^M3) 

1000 \ 




100 200 300 400 500 
Ml (GeV) 



Figure 6: Present upper bound on Xm(e~*'^*''A/i3). The reference values tan/? = 10, 
A/Mr = 10^ and A/M3 = 2 ■ 10'^ have been taken. Solid (dashed) lines correspond 
to case a) (b)). 

Let us firstly discuss the range. For our reference values, the upper estimate 
for Im{6^^S^^)ii is O(10~^) and the present bound on de already constrains it to 
be smaller in a large region of the plane. Although we already know that they are 
not dominant, it is worth to discuss the LR RR amplitude and the similar LL LR 
one because, as already mentioned, they are quite model independent. Allowing for 
higher triplet masses, however, these amplitudes shift below the present experimen- 
tal sensitivity. This is not the case for the LL RR amplitude. The maximum value 
allowed for 2m{6^^mi6^^)ii/mr is displayed in the right panel of fig. El for case a). 
In case b) it is ~ 10~^ everywhere. Then, the present experimental bound on de has 
already explored the radiative range for roughly 2 — 3 orders of magnitude. 

This can be translated into an upper limit on Xm(e~*''^"'A/i3), as in fig. IHl Notice 
that this limit is indeed very strong and can be hardly evaded even looking for less 
favorable parameters than those taken as reference. To check this, keep e.g. 
fixed and try to worsen the limit: a reduction of /i with respect to is unlikely 
to reduce it more than one order of magnitude; a suppression by another factor 10 
would require a A/Mt ~ 1.6 ; on the other hand, for values of tan/5 larger than 10 
the limit linearly improves. As a result, in the framework of the seesaw accompanied 
by SU{5), the limit on Jm(e-*'^"A/'i3) can be considered robust. 
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Experiment 
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Figure 7: Same as fig. El but for < 10^^*^ e cm. 



Theory: SU (5) with Seesaw 
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Figure 8: Same as fig. El but for i = 2. 
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Upper limit on 
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Figure 9: The upper bound on Xm(A/23) which would be extracted by improving the 
present limit by many orders of magnitude ^l^j. Solid (dashed) lines correspond 
to case a) (b)). 



Let us now turn to discuss d^. By comparing figs. [Z|and|Hl it turns out that the 
upper bound on the LR RR insertion would require a sensitivity to at O(10~^^) 
e cm, at the very limit of planned experimental improvements. Quantitatively, the 
upper estimate for Im{5^^ 771^5^^)22 /tt^t which, for case b) is ~ 4 10~^, is more 
promising. The major part of the plane in fig. |H1 could be tested with at the 
level of 10~^^-10~^^ e cm and in general ci^ should not exceed O(10~^'^) e cm. As 
a result, the eventual presence of triplets doesn't enhance by much the range for 
d^ with respect to the pure seesaw case. Nevertheless, here too the possibility of 
constraining Xm(A/23) can be envisaged, as shown in fig. IHl Notice that, due to (1221), 
a limit on Xm{M2z) comparable to the present one on Xm(e~*'^"A/i3) would require 
to improve the d^ sensitivity down to 5 10^^'^ e cm. 
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5 Conclusions 



Planned experiments might significantly strengthen the limit on f21 ^] and 
^lE]. Their eventual discovery could be interpreted as an indirect manifestation 
of supersymmetry but could not reveal which source of CP (and possibly flavour) 
violation is actually responsible for the measured effect. Clearly, all sources in 
principle able to give the lepton EDM even at an higher level would be automatically 
constrained while those which fail in giving the lepton EDM at the desired level 
would be automatically excluded from the list of possible candidates. 

In this work, we have estimated the ranges for the lepton EDMs induced by the 
Yukawa interactions of the heavy neutrinos, both alone and with the simultaneous 
presence of the heavy SU (5) triplets. It turns out that the FV LLRR amplitude is 
in general larger or comparable to the PC one. So, EDMs are enhanced for large 
values of tan/5 and do not strongly depend on oq. 

The pure seesaw, even with large tan (3 and very hierarchical right-handed neu- 
trinos, cannot account for d^ above 10"^^ e cm. Its eventual discovery above this 
level would then signal the presence of some source of CP and LP violation other 
than the neutrino Yukawa couplings. The heavy triplets Yukawa couplings would be 
excluded from the list of possible sources because their additional presence do not 
significantly enhance the predicted range for d^. Notice however that in the latter 
case a hierarchical right-handed neutrino spectrum is no more essential to end up 
with d^ at an interesting level for planned searches. Prom the theoretical point of 
view, finding li^ above 10~^^ e cm would indeed have a remarkable impact. 

Interestingly enough, the present experimental sensitivity to d^ is already testing 
the simultaneous presence of triplets and right-handed neutrinos. Correspondingly, 
constrains on Xm(e~*''^"A/i3) have been derived which are significant even for quite 
large values of the triplet mass and unfavorable supersymmetric masses. Without 
the triplets, the radiatively-induced de is close to the present experimental sensitivity 
only in models with large tan (3 and small slepton masses. Therefore, an experimental 
improvement would eventually provide interesting limits on the imaginary part of 
the relevant combination of neutrino Yukawa couplings and right-handed neutrino 
masses. 

In the present discussion we have been looking for results as general as possible. 
Indeed, the specification of any particular seesaw model has been avoided and the 
attention has rather focused on the dependencies on the supersymmetric masses 
and heavy thresholds. Although some relevant seesaw models deserve a dedicated 
analysis j2H], figs. EUSlandlHlare suitable for a quick check of the status of any given 
seesaw model with respect to the present and planned experimental limits on lepton 
EDMs. 
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A Appendix: Seesaw 



In the basis where charged fermion and right handed Majorana neutrino masses are diagonal 

W 3 u^^yuQHu + (f^VdQHd + e^^yeLHd + v^^yuLHu + Mrv"^ (33) 

where Q = {u d)'^ , L ~ {v e)^ and {H^f^^yi = Vd{u)- Soft scalar masses are defined as 

Csoft 9 ulfinluR + d\^mldB, + Q^'mQQ + e\mleR + Vm\L + v'^^m'lvR 

+ + JJ^AddiWd + eJ^^eeiUd + i'jj^i.i'LWii + /i-c.) (34) 

Let us also introduce the following notations: 

ylVx = X y^yl = X {x = e,u, d) 

Pay. = y'~^'> ylPay. = Na Pay.ylPa = Na (a = 1,2, 3) 

where P2, Pi project out M3 and M3.2 respectively, P2 ~ diag(l,l,0). Pi = diag(l,0, 0), and 
Pz = i3. 

A.l Running 

Defining t = j^pp InQ, the running of the Yukawa coupling constants is governed by: 

^ = yi'^^iNa + E + D^^)] 



dye 
dt 

dyu 
dt 

dVd 
dt 



ye[iE + Na + De] 

yu[W + D + D'^^^] (35) 
yd[iD + U + Dd] 



dt 
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where 



= [Tr(3C/ + 7V„)-(35i + ?.9?)]i3 

o 

= [Tr(3D + i?) - (85! + |g?)]i3 

= [Tr{3U + N,)-{f9i + 3gl + ^gl)]i3 (36) 
D, = [Tr{3D + E)-{fgl + 3gl + ^gl)]i3. 



For the trilincar couphngs, defining PaAi, = A^f^ 



dAe 

dt 
dAu 

dt 

dAd 
dt 



4^Ae + 5Ae^ + 2y^y^^^^ Ai"-^ + A^N„, + D,A, + 2D,y, (37) 
4UAu + 5A^U + 2yuy\Ai + A^D + D^'') A„ + 2b<^^y^ 
AD Ad + f>AdD + 2yayiAu + A^U + DaAa + 2bdyd 



where 



^(«) = [Tr{3ylA^ + ^/l^^^l")) - {3glM2 + plMr)V3 
b, = [Tr{3y\Ad + ylA,)-{3glM2 + \glm)V3 (38) 

^(«) = [Tr(32/iA„+y('^)U("))-(^532M3 + 3ff2M2 + H52^^)]^^ 

ba = [Tr(3t/>d + ylA,) - (^^1^3 + 3ff|M2 + ^5?Mi)]i3 . 

For soft scalars, defining Pam^Pa = m1^°'^ 

^ = {mi, E + N^] + 2{ylmlye + m|,^£ + AIA,) + 2(y(<')^m2("^t/(«) + 7V„ + A^-)^ A^^'^) + Gl 
dml 



dt 
dm^ 
dt 

dt 

dt 

dt 

where 



2{ml E} + Hyemhl + m%^E + A,Al) + Ge 
2{m^^"\ 7VJ + 4(y(«)mi2/W^ + ml,„iV„ + A^AI")^) (39) 
{ml, U + D} + 2{yimlyu + mljU + AIAu) + 2(y^m^yd + m]j^D + A^Aa) + Gq 
2{ml, U} + 4(t/„m|2/t + m%^U + A^Al) + G„ 
2{ml b] + 4(t/dm|yJ + m%^b + A^A^) + G^ 



Gl = -ilglM? + 6glM^)i3 G, = -{'^glMl)i3 

2 32 
Gq = -(Y^5?Mf + GglMl + -glMi)i3 (40) 

G„ = -i^glM? + f 53^^3^)13 Gd = -i^glM? + f fflM|)i3 - 
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dt 



Finally, 

= QTr{ydmly\ + y\mlyd + ml^D + A^Ad) + 2Tr{y,mlyl + ylmly^ + m]j^E + AIA,) + Gl 
= 6Tr(2/„m|yt +yt^2y^ + ^2^^t/ + ^tA„) (41) 

and we define m^'rr = rfm?r /c?i. 

For energy scales between A and -M3, one has to take a — below Afs and above M2, a — 2] 
while below M2 and above Mi, a = 1. 



B Appendix: SU{b) + See-saw 

We adopt the following notation to write matter and Higgs superfields: 



H3d2 

\-H2fJ 

where < H^d(2u) '"d{u)- The superpotential 
1 



/ 





u§ 


-U2 








/ 


dl \ 

dl 









ul 




d2 1 








"2 


~ui 





vP' 




05 = 




dl 



















e 


V 


-d^ 


-d^ 


-d^ 




0^ 









(42) 



H 



TT 2 \ 

H2u^ 
\h2u"J 



(43) 



1 



W3- ^""""yu ^''"H^'tABCDE + V2 i^^^y, 4>aHb + my. <I>aH^ + 2^1^^ ^ (44) 

where A,B, ... — 1, ...,5 and flavour indices are understood, gives rise to (|33|l with y^ — y^ and 
Uu = Uu- For the soft breaking part of the Lagrangian 



+ ( u'^'^AuU Vu + d"^ Ai d Vd + e"-" A^e Vd + u"^ A^v u„ ) + h.c 



JcT aT 



(45) 



where the scalar fields are ^ = {u'^ , u, d, ff^) , cf) = {d'^^e^v) and f) = u'^ and gauginos are denoted 
with A5. In this way, in the scalar lepton mass matrix, ra\j^ = to^, m^jui = ra^ . 

B.l Running 

Setting t EE InQ, dg^/dt = -Syf, dM^/dt = -63^/5 and 



dt 
dyu 
dt 



6yeE + Wye + yeN + G^ye 

uVu 

6y^N + Ay^E + G 



(46) 
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where = -84/5g| + ATr{E), G„ = -96/5.g| + Tr{W + N),G^^ -48/55^ + Tr{W + N). For 
scalar masses 

= {ml*,2E + 3U} + ^yemlyl + mlE + A,Al) + 6{y^mlyi + m,lU + A.,Ai) + G^ 
= {ml , AE + N} + 8{ylml*y, + m|Z? + 4^) + 2{ylml*y, + m^iV + ^A,) + G^ 

dm?* 

= b{ml\N} + my,mlyl + mlN^A,Al) (47) 

where G^, = -144/5c/|M|i3, = -96/5g|M|i3. For triUnear couphngs 
dA 

^ 10EA, + WA, + 8A,E + &A^yiy, + A,N + 2y,ylA, + G,A, + 2Geye (48) 

^ = 9?7A„ + 2^A„ + 9A„[/ + 4Ae2/ly„ + 2Ai? + 4y„yt^<i + G„A, + 2G„y„ (49) 
dA 

= 7A^A, + llA,7V + 4A,£; + 82;,2/lv4e + G,A, + 2G,y, (50) 

where G,. = -84/5ffiM5 + 4Tr(ytAe), G, - -96/55|M5 + Tr(3j/tA„ + zjt^,), G, = -48/55pf5 + 
Tr{iy\Au + y^Ay). Finally, for scalar higgses 



dt 



Tr{QU + 2N)ml + QTr{y^mlyl + ylm^^yu + A«4) 

96 

+ 2Triy,mlyl + y^my^ + A,Al)--glMl (51) 
^ = 8Tr{E)ml + 8Tr{y,mlyl + y^^y, + A,A\) - ^.g^Af ^ . 
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